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Abstract. In a previous paper [1] we solved a countably infinite family of one- 
dimensional Schrodinger equations by showing that they were supersymmetric partner 
potentials of the standard quantum harmonic oscillator. In this work we extend these 
results to find the complete set of real partner potentials of the harmonic oscillator, 
showing that these depend upon two continuous parameters. Their spectra are identical 
to that of the harmonic oscillator, except that the ground state energy becomes a 
tunable parameter. We finally use these potentials to analyse the physical problem of 
Bose-Einstein condensation in an atomic gas trapped in a dimple potential. 
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1. Introduction 

In a previous paper in this journal [1], we found the exact solution of the Schrodinger 
equation for a family of potentials consisting of a harmonic oscillator term plus a rational 
function. The Schrodinger equation for these potentials is 



+ 



%{xf-Hp{x)H';{x) 



where the 'Hp{x) are the pseudo-Hermite polynomials defined by 



Hpix) 



-iyHp{ix) 



dxP 



These pseudo-Hermite polynomials, 'Hp{x), are just the regular Hermite polynomials, 
Hp{x), with all the minus signs removed. For even p the potentials take the form of 
a harmonic well with an additional attractive dimple of depth Ap at the origin. The 
bound states have energies, = 2k + 3, and wavefunctions. 
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Hk+i{x) + 



'Hp{x 



-Hkix) 



2^+m{k + p + i)y/^ 

where = 0, 1, 2, 3, ... In addition, there is a new ground state with energy, E_i = 1- 
and wavef unction, 
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It is as if the attractive dimple of depth Ap has taken the ground state energy level of 
the harmonic oscillator and lowered it by an amount 2p, whilst leaving all the other 
energy levels fixed. In the above solution, p has to be an even integer, but one could 
imagine a situation where p is an arbitrary positive real number. This would yield a 
potential with a spectrum consisting of an evenly spaced ladder with a ground state 
energy level which may be lowered at will by the addition of an attractive dimple. 
Such a potential, as well as being interesting from a purely mathematical point of view, 
may have practical application in the area of atomic Bose-Einstein condensation. It 
has been established, both theoretically [2] and experimentally [3], that the adiabatic 
addition of a narrow attractive potential to a harmonic trap can lead to Bose-Einstein 
condensation at a higher critical temperature. The theoretical justification for this 
effect is that an additional bound state may be formed in the superposed potential, 
which may have energy as low as the chemical potential of the unperturbed system, 
leading to macroscopic occupation of that level. The problem with such a heuristic 
explanation is that when two potentials are superposed, one does not usually obtain a 
spectrum which bears any resemblance to that of either of the superposed potentials. 
Our potential is therefore special in that the effect of adding an attractive dimple of the 
appropriate form has the sole effect of lowering its ground state energy. This makes it a 
natural starting point for future theoretical work on the dimple problem in atom traps. 

In this paper we generalize our previous family of soluble potentials to allow the 
parameter p to take non-integer values, and thus obtain an energy spectrum with a 
continuously tunable ground state. Since our previous potentials were solved by noting 
that they are supersymmetric partner potentials of the harmonic oscillator, we proceed 
by determining the most general potential which is a partner of the harmonic oscillator. 
This leads to a two-parameter family of soluble potentials, of which the desired potentials 
are a subset. 

The rest of the paper is organized as follows. In section 2 we give a prescription for 
generating a two-parameter family of soluble potentials from a known soluble potential 
using the factorization approach. In section 3 we apply this prescription to obtain 
the set of all potentials which are partners of the harmonic oscillator, and derive the 
eigenvalues and eigenfunctions for these potentials. In section 4 we plot a representative 
sample of these potentials and their eigenfunctions. In section 5 we analyze the problem 
of Bose-Einstein condensation in one of the dimple potentials derived in section 3, and 
in section 6 we make our summary and conclusions. 

2. The factorization approach 

Consider the one-dimensional single-particle Schrodinger equation. 




A New Two-Parameter Family of Potentials with a Tunable Ground State 



3 



The idea of the factorization approach (see |1] and references therein) is to write the 
Hamiltonian operator, Hi, as a product of two first-order differential operators, 

A^A, 



where 



so that 



Hi 



A 



d_ 

dx 



+ W(x), A^ = -^ + W(x), 

dx 



Vi{x) = W\x) - W'{x). 

If we now define the operator H2 = AA'^ by reversing the order of A and A\ we see that 
H2 is a Hamiltonian operator corresponding to a new potential, V2{x), 

d^ 



Ho = AA^ 



dx^ 



+ V2{X) 



V2{x) = W{xy + W'{x). 



The potentials Vi{x) and V2{x) are known as supersymmetric partner potentials. 

The eigenvalues and eigenfunctions of Hi and H2 are related. If ip^^ is an 
eigenfunction of Hi with energy eigenvalue El^\ then Aip^^ is an eigenfunction of H2 
with energy eigenvalue E^'' since 



Ho 



A 



Similarly if if)^'' is an eigenfunction of H2 with eigenvalue then A'^t/j^^'^ is an 

eigenfunction of Hi with energy eigenvalue E^^^^ since 



Hi [aV. 



(2) 



At 



H2^ 



(2) 



(2) 



It follows that if we know how to exactly solve one of Hi or H2, we can immediately 
derive the exact solution of the other. 

Suppose that we know the exact solution of the Hamiltonian H2. Let us find its 
complete set of partner Hamiltonians, Hi. The function W{x) in the factorization must 
satisfy the equation 

W{xf + W'{x) = V2{x)-X, 

where we have noted that the addition of an arbitrary constant A does not affect the 
solution of H2. To solve this equation we make the standard substitution 

iy(x) = -^ln0(a;)-^ 
dx 



which yields the result 
_d^^ 
dx 



(f){x) 

2 +V2{x)(j){x) = X(f){x). 



It follows that (j){x) is a solution of the original Schrodinger equation for the known 
soluble potential V2{x). The central idea of our approach is that whilst a 

solution of the Schrodinger differential equation, does not have to be a normalizable 
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wavefunction. There is therefore a general solution, (f)x{x), for any value of the pseudo- 
energy parameter, A, 

where (l)x,i{x) and (f)x,2{x) are two independent solutions of the Schrodinger equation. If 
V2{x) is even, these solutions will have definite parity, and without loss of generality we 
can assume (px,i{x) is even and 0A,2(a^) is odd. The function W{x) is then given by 

0A,lW + S(px,2{x) 

where s = D/C. We see that Wx,s{,x) depends upon the two real parameters A and s, 
and therefore so does Vi{x). The only restriction we will impose on A and s is that the 
function (l)x{x) has no zeros, so that the potential Vi{x) is non-singular. This means 
that A must be less than or equal to the ground state energy of H2. 



3. The complete set of partner potentials of the harmonic oscillator 

In this section we will take the harmonic potential, V2{x) = a;^, to be the known soluble 
potential, and derive the complete set of partner potentials, Vi{x). 

3.1. Schrodinger equation for harmonic oscillator 

To solve the Schrodinger equation for the harmonic oscillator. 



+ X (j) = X(f), 



dx^ 

we first remove the asymptotic behaviour by setting 

0(x) = e(x)e--'/^ 

which yields the equation 

^|-24 + (A-l){ = 0. 

If we change the independent variable to z = x"^, the latter equation becomes, 

which has the form of the confluent hypergeometric equation [5], 

z— + {c-z)—-a^ = 0, 
dz"^ dz 

with a = |(1 — A) and c = |. This can be solved using the method of Frobenius to give 
the two solutions 

$(a, c;z) and ^-^""^$(0 + 1 - c, 2 - c; z), 

where 

u ^ ^ r(a + n)r(c)^" 
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It follows that the two independent solutions of the harmonic oscillator Schrodinger 
equation are 

To obtain the standard harmonic oscillator eigenfunctions, tlj^^\x), we write X = 2k + 1 
and note that the confluent hypergeometric function, $(a, c; z), becomes a terminating 
polynomial if parameter a is a negative integer. If k is even, then the even solution 
terminates to give 

whilst if k is odd, then the odd solution terminates to give 

These can be written in terms of Hermite polynomials using the identities j5], 

H2n{x) =(-l)"M:$(_n,|;x2) 
nl ^ 

so that the normalized solutions of the harmonic oscillator are 

1/2 

In our previous paper, we found suitable nodeless solutions, by employing the 
internal symmetry x — )■ ix, A — > —A in the harmonic oscillator Schrodinger equation. 
Our two general solutions then become 

if p is an even or odd integer then (/>p,i(x) or 0p 2(a;), respectively, are the solutions 
employed in our previous paper. The point of what we are doing here is to ensure 
that our notation aligns with that of our previous paper for integer p; however we will 
henceforth be considering the general case where "p may be non-integer. Note that the 
confluent hypergeometric identity [5], 

$(a, c; —z) = e~^$(c — a, c; z), 

implies that 

= $(2±1, i; a;2)e-^'/^ <P,,2{x) = x$(f + 1, i; x')e-^'/', 
which we recognize as the general solutions for A = — (2p + 1), as we would expect. 



(2), 



X) 



2^k\^ 
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3.2. Even parity partner potentials 

To construct new soluble potentials, Vi{x), we first evaluate the superpotential W{x). 
Let us first restrict ourselves to the even potentials by using the solution (pp,i{x) so that 

d 



W{x) 



dx 



ln0p^i(x) 
d 



ln$(£±i,l; 

dx 

—X + 2{p + l)a; 



X 



2 ' 2' • 



2 ' 2' ■ 



2 ' 2' ■ 



-$(a + l,c+ 1,2;). 



where we used the identity [5], 

From this we generate 

V2{x) = W^{x) + W\x) = x^ + 2p + I = Vnix) + 2p + 1, 

which is clearly the harmonic oscillator potential, Vh{x), with the offset 2p+ 1. We will 
choose the offset for Vi{x) so that it asymptotically approaches Vh{x) at infinity. This 
leads to an offset of 2p — 1. Hence 

Vi{x) = W\x) - W'{x) = 2W\x) - Vi{x) = Vp{x) + 2p-l, 

giving the final result for the partner potential 



Vpix) 



x^ — Axw{x) + 2w'^{x) — 4p, 



where we define for convenience 



w[x] 



2{p+l)x 



V 2 ' 2 ' ' 



The normalized eigenf unctions of this potentials are then given by 

1 ni/2 



(2) 
k J 



(2), 



X] 



n V2 r 



1 



X + w[x) 

dx 



HJx)e-^"'^ 



1/2 



w(x)Hk(x) -2kHk-i(x) 



-x2/2 



_2''+^k\{k + p+l)y^_ 
with energy E^ = 2k + 3 for k = 0, 1, 2, 3, . . .. The new ground state is given by 

1 p^V2 

ip^i{x) oc 



$(-2+1 i-a;2V 
V 2 ' 2 ' ' 



The normalization factor can be guessed by analytically continuing from the result for 
even integer p derived in the previous paper [1] , 



p\ 2P 



1/2 -a;2/2 



Hp{x) 
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and using the result that 



1-Lp{x) 



\ 2 ' 2 ' ■ 



This gives the final result for the normalized ground state wavefunction, 

2P 



r(p+l)v^. 



2 ' 2' ' 



with energy, E^^ = 1 — 2p. We see that the spectrum of Vp{x) is identical to that of 
the harmonic oscillator except that the ground state energy has been lowered by the 
amount 2p. 

Notice that in the above we are guessing the normalization integral, 



V 2 ' 2 ' 



dx 



2Pr(2±i)2 ' 



ip > -1) 



which we previously proved pLj for even integer p. A formal proof of this result is given 
in Appendix A. 



3.3. General (skewed) partner potentials 

Let us now consider the most general non-singular partner potentials, Vi{x). We do this 
by using the general solution. 



X) + S(l)po(x) 



V 2 ' 2 



2' ' 



-x'2/2 



where the skewness parameter, s, determines how much odd solution we add to the even 
solution. We choose s so that 4>p{x) has no zeros, and hence Vi{x) has no singularities. 
To determine the maximum allowed range of s, we use the result for the asymptotic 
limit of the confiuent hypergeometric function [5], 



^(a, c, z) ~ 



Tic) 



r(c — a) 



T(c) 



r(a) 



\z\ — )■ 00 



so that 



■'p,s 



X] 



r(2±i) r(f + i 



-sgn(xj 



In order for to not change sign, we require the condition 

_ ] 

\s\ *\ 5 rn.n.T. 2 " 



r (f + i) 



Note that we need only consider < s < Smaxi since the potentials for negative s are 
just the mirror images of those for positive s. All the formulae for the eigenvalues and 
eigenf unctions are the same as in the even case, except that w{x) is now given by 



w[x] 



2(p+ l)x$( 



0+3 3. ^2 

2 ' 2' 



x^) + sx^il + 1, f ; x^) + \{p + 2)sx2$(f + 2 



5. ^2 



' 2' 



X 



<l>(2±i,i;a;2) + sx$(| + l,f;x2) 
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The new ground state is given by 



To normahze this we need to evaluate the normahzation integral 



-c 



$(2±i i;a;2) + sx$(f + 1 



3. ^2 



' 2' 



■ dx. 



We can guess the value of this integral by noting that we know its value at s = 0, and 
it must diverge at s = is^ax, to obtain a sensible interpolation formula 

ip,o _ xAFr(p + i) 



Ip,s — 



2P 



r(| + i)2-fr(£±i)2 



A formal proof of this result is given in Appendix A. 



3.4- Summary 

Wc have derived the most general non-singular partner potentials of the harmonic 
oscillator in terms of two continuous real parameters p and s. and have presented their 
eigenvalues and eigenfunctions. All these potentials have the same spectrum as the 
harmonic oscillator, except that the ground state has been lowered by an amount 2p. 
The parameter s is a skewness parameter which changes the shape of the potential 
without changing its spectrum. The allowed ranges of the parameters p and s are 



P ^ 1) I'^l ^max — 2 



r(f + 1) 



r(2±i) • 

Note that for — 1 < p < 0, the ground state energy is raised and can be brought 
arbitrarily close to the first excited energy level; we find in this case that the potential 
has a positive bump at the centre, and so is a double well. In the following section we 
will plot the potentials and their eigenfunctions for representative ranges of p and s. 



4. Representative examples of our new soluble potentials 

In the previous section we derived the complete solution of the Schrodinger equation 
for a two-parameter family of potentials. We now wish to plot the potentials and their 
eigenfunctions for suitable regions of the (p, s) parameter space. For ease of notation 
let us henceforth rescale s — >■ s/s^ax- 



4.1. Symmetric dimple potentials: p > 0, s = 

If the skew parameter s equals zero, the partner potential generated will be even. If 
p > 0, the potential will take the form of a harmonic well plus an attractive dimple of 
depth 4p at the origin. The potentials and their first three wavefunctions are shown 
in Fig. (1) for p = 0.5, 1, 2 and 4. We see that the ground state becomes narrower 
and taller as the dimple depth increases, whilst the excited states are less affected. It 




X 



Figure 1. Symmetric dimple potentials and their wavefunctions corresponding to 
parameter values p > and s = 0: (a) Potentials; (b) Ground state; (c) First excited 
state; (d) Second excited states. 

is as if only the ground state really notices the attractive dimple, and has its energy 
pulled down by half the depth of the dimple. The other states keep the same energy 
and change their shape only enough to retain orthogonality to the ground state. 

4-2. Symmetric double well potentials: > p > —1 , s = 

If > p > — 1 the potential takes the form of a double well sitting at the bottom of a 
harmonic potential. The potentials and their first three wavefunctions are shown in Fig. 
(2) for p = —0.3, —0.5, —0.7 and -0.9. The double well shape can partly be understood 
as the effect of a "negative dimple" which pushes up the centre of the potential to form 
a barrier between two wells. As p approaches —1 the centre of the potential becomes 
flatter and then develops a minimum which resembles the original parabolic minimum, 
but raised up by 4 units of energy, whilst the two wells move up the sides of the harmonic 
potential. The ground state starts off with a single maximum which then flattens out 
and becomes a minimum between two maxima. The wavefunction initially straddles 
the two wells when they are relatively shallow, but as they become deeper and further 



Figure 2. Symmetric double well potentials and their wavefunctions corresponding 
to parameter values > p > — 1 and s = 0: (a) Potentials; (b) Ground states; (c) First 
excited states; (d) Second excited states. 

apart, the wavef unction splits into two pieces centred on the two wells. The first excited 
state wavefunction also initially straddles both wells when they are shallow, and splits 
into two pieces as they become deeper. In the extreme limit where p approaches — 1, 
the ground state and first excited state take the form of symmetric and antisymmetric 
combinations of two peaks centred on the two wells. This is to be expected since the 
lowest two energy levels are now very close together. The second and higher excited 
states are relatively unaffected by the double well potential, and their wavefunction 
always straddle both wells. The lowest two energy levels effectively form a two-level 
system with a large energy gap to the remaining harmonic ladder. 

4-3. Skewed dimple potentials: p > 0, 1 > s > —1 

If p > and 1 > s > 0, the potential takes the form of an attractive dimple whose 
minimum is skewed to the right (for > s > — 1 the mirror image is obtained). The 
potentials and their first three wavefunctions are shown in Fig. (3) for p = 4 and 
s = 0, 0.9, 0.99 and 0.999. The amount by which the dimple moves to the right depends 
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Figure 3. Skewed dimple potentials and their wavefunctions corresponding to 
parameter values p = 4.0 and 1 > s > 0: (a) Potentials; (b) Ground states; (c) 
First excited states; (d) Second excited states. 

logarithmically on 1 — s. The ground state moves with the minimum of the dimple, 
and roughly maintains its gaussian shape. The other excited states do not move with 
the dimple, and change their shape to retain orthogonality to the ground state. The 
first excited state starts off with odd parity when s = 0, and tends asymptotically to 
even parity as s approaches 1 and the dimple moves off to infinity. This is because 
the orthogonality to the ground state becomes less restrictive as the ground state 
wavefunction moves off to infinity with the dimple. The second excited state similarly 
evolves from even parity at s = to odd parity as s approaches 1. 

4 ■4- Skewed double well potentials: > p > —1, l>s>— 1 

If > p > — 1 and 1 > s > 0, the potential takes the form of a double well with the 
right hand well deeper and further from the origin than the left (for > s > — 1 the 
mirror image is obtained) . The potentials and their first three wavefunctions are shown 
in Fig. (4) for p = —0.9 and s = 0, 0.9, 0.99 and 0.999. The amount by which the deeper 
well moves to the right, and its depth, depend logarithmically on 1 — s. In contrast the 





Figure 4. Skewed double well potentials and their wavefunctions corresponding to 
parameter values p = —0.9 and 1 > s > 0: (a) Potentials; (b) Ground states; (c) First 
excited states; (d) Second excited states. 



shallow well rapidly approaches a limiting form. The ground state wavefunction has 
two peaks which are centred on the two wells. As s increases, the height of the peak 
in the right hand well increases, whilst that in the left decreases. As s approaches 1, 
the ground state looks like a gaussian centred on the right well. The first excited state 
starts off with odd parity at s = 0, with equal sized peaks in two wells. As s increases 
to 1, the peak in the right hand well vanishes, and the wavefunction looks like a skewed 
gaussian centred on the left hand well. It is as if the right hand well has moved so far 
away that the first excited state is effectively the ground state of the left well. The 
second excited state starts off with even parity at s = 0, and ends up looking like the 
first excited state of the left well as s approaches 1. 

5. BEC condensation temperature in a dimple trap 

Although Bose-Einstein condensation is not possible in one or two dimensions in a 
translationally invariant system, no such restriction exists in a confined system [6]. 
Low-dimensional condensates are routinely produced in magnetic or optical traps that 
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are well modelled by a harmonic oscillator potential. The critical temperature, T^, for a 
cloud of bosons in such a one-dimensional harmonic well can be found from particle 
number conservation when the chemical potential equals the ground state energy. This 
gives the equation 

where /3c = l/ksTc- For large this approximates to 

N 



^BJ-c 



■ hu. 



InN 

Experimentally the condensate fraction may be increased by adiabatically deforming 
the shape of the trapping potential [3] , so that it takes the form of a harmonic well with 
an attractive dimple at the origin. Such dimple potentials are highly reminiscent of the 
potentials discussed in section 4.1. These potentials allow condensation to be achieved 
at a higher temperature, and in a reversible manner. Although heuristic semiclassical 
explanations of this phenomena exist in the literature [2], it is preferable to have an 
exactly soluble quantum model. We will therefore consider Bose-Einstein condensation 
in a dimple trap of the form discussed in section 4.1, with depth 4p{huj/2). The energy 
levels of this potential are Eq = {1 — 2p)huj/2 and En = {2n + l)huj/2 {n = 1, 2, 3 . . .), 
so that the critical temperature is given by 



n=l 



The dependence of Tc as a function of p is plotted in Fig. (5) for a fixed number of 
atoms, A^ = 10^. As expected, increases monotonically with and doubles when 
p ^ 200. This corresponds to a dimple depth of Ap{h,uj/2) = iOOhu. This is easily 
achievable in experiments [7], where typically w ~ 27r x lOHz so that AOOhu ~ 2nK, 
whilst dimple depths of up to lOOOni^ may be produced. 



6. Summary and Conclusions 

In this paper we have derived the most general non-singular partner potentials of the 
quantum harmonic oscillator, and obtained their complete set of energy eigenvalues and 
normalized eigenfunctions. These potentials depend upon two continuous parameters: 
a depth parameter, p, which tunes the depth of the centre of the potential, and a skew 
parameter, s, which breaks the reflection symmetry of the potential. The spectra of 
these potentials are identical to that of the quantum harmonic oscillator, except for the 
ground state whose energy is lowered by 2p. A possible physical application of these 
potentials is to the investigation of Bose-Einstein condensation in dimple traps. 
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Figure 5. Relative transition temperature, Tc{p)/Tc{0), as a function of dimple depth 
parameter, p, for a cloud of N = 10^ bosons in a dimple trap. 



Appendix A. Normalization of ground state eigenfunctions 



Suppose that (j)p^si{x) and (j)p^s2{^) independent solutions of the Schrodinger 
differential equation, 



+ V{x)(f) = epcf). 



It follows that the Wronskian, 

W{x) = (t>p,sA^)(t>'p,s2i^) - (f)p^,^{x)(j)p,s2{x), 

is a constant, which we will denote by C(p, Si,S2). Dividing the Wronskian by 0^ 
gives 



SI 



x 



dx y(j)p,s-,{x) 



4>p,sAxY (f>p,sAx) (t>p,sAxY 

so that upon integration we get 

dx _ 1 [" 0p,^.,(.t) 

For our problem 

(t>p,s{x) = [$(£±i, i; x^) + sx^il + 1, 1; x^) 
which has asymptotic limit as x — >■ ±oo 

r r(i) r(^) 

Evaluating the Wronskian at x = gives C(p, Si,S2) = S2 — Si, so that the integral 
becomes 



°° dx 



1 



oo <pp,sr{xy S2-S1 



a + /3s2 a — /3s2 
a + a — (3si 



2a(5 
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where for convenience we have defined a = r(|)/r(2±i) and /3 = r(|)/r(| + l). It 
foUows that 

dx r(| + i)r(^^) yir(p+i) 

„2, 



. <l>p,s{^r r(| + 1)2 - f r(£±i)2 2P [r(| + 1)^ - f r(2±i)2] ' 

where in the last step we used the duphcation formula 2'^^^^r{z)r{z + ^) — \^r{2z) 
with z — (p + l)/2. We see that the simple formula obtained by guesswork based 
on analytic continuation from even integer p, and the divergence of the integral at 
s = ±Smax, is correct. Moreover, the Wronskian technique demonstrated above allows 
one to normalize any new ground state wavefunction generated by the factorization 
method. 



References 

[1] Fellows J M and Smith R A 2009, J. Phys. A: Math. Theor. 42, 335303. 
[2] Pitaevskii L and Stringari S 2003, Bose-Einstein Condensation (Oxford: OUP) pp 158-160. 
[3] Stamper-Kurn D M, Miesner H-J, Chikkatur A P, Inuoye S, Stenger J and Ketterle W 1998, Phys. 
Rev. Lett. 81, 2194. 

[4] Cooper F, Khare A and Sukhatme U 2001, Supersymmetry in Quantum Mechanics (Singapore: 
World Scientific). 

[5] Gradshteyn I S and Ryzhik I M 1980, Table of Integrals, Series and Products (London: Academic). 
[6] Petrov D S, Gangardt D M and Shlyapnikov G V 2004, J. Physique IV France 116, 5. 
[7] Garrett M C et al 2011, Phys. Rev. A 83 013630. 



